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A set of novel derivative terms for spin-2 fields were proposed. They are the wedge products of curvature 
two-forms and vielbeins. In this work, we investigate the properties of novel two-derivative terms in the context 
of bi-gravity. Based on a minisuperspace analysis, we identify a large class of bi-gravity models where the 
Boulware-Deser ghost could be absent. We give a new perspective that Weyl Gravity and New Massive Gravity 
belong to this class of bi-gravity models involving novel derivative terms. In addition, we discuss the UV cut-off 
scales, dynamical symmetric conditions and novel higher-derivative terms. 


I. INTRODUCTION. 

In the pursuit of going beyond Einstein, extensive efforts 
have been devoted to constructing gravitational theories that 
are different from general relativity. It is believed that the 
Einstein-Hilbert action with a cosmological constant is the 
only consistent nonlinear action for single massless spin-two 
held, so additional ingredients are necessary: Eierz-Pauli the¬ 
ory gives the graviton a mass term lH; Brans-Dicke theory 
introduces a scalar degree of freedom ||2l ; Lovelock theory 
involves higher spacetime dimensions and higher derivative 
terms Jl]. 

Motivated by the unexpected accelerating expansion of 
the present universe |01, more models were constructed in 
recent attempts, such as high derivative scalar theories iH] and 
nonlinear massive gravity i^. One of the guiding principles 
is that a consistent model should be free of Ostrogradsky’s 
ghost arising from higher order equations of motion. 

Antisymmetrization is a universal element in these new 
models. Based on this ingredient, a general framework 
was developed for ghost-freeQ Lorentz-invariant, Lagrangian 
held theories ||3,lltl - In this framework, a set of novel kinetic 
terms for spin-2 helds were proposed 

where d is the number of spacetime dimensions. They are the 
wedge products of geometric differential forms: a curvature 
two-form and several vielbeins. The vielbeins can be the 
same or different^. Geometric intuition was used to construct 
these nonlinear terms. 


* lii.wenliang@gmail.com 

^ Let us clarify that "ghost-free" in this framework means the building blocks 
are potentially free of Ostrogradsky’s scalar ghost. In other words, the 
con'esponding equations of motion for the scalar modes could be at most 
of second order. The scalar modes may come from the Helmholtz-Hodge 
decomposition of tensor fields. 

^ If all the vielbeins coincide, we have the standard Einstein-Hilbert kinetic 
term. 


Along the line of Lovelock theory, we can build novel 
higher-derivative terms for spin-2 helds by introducing more 
curvature two-forms into the wedge products 

A • • • A A E^^'> A • • • A E^'^\ (2) 

which is possible when spacetime has more than four dimen¬ 
sions. Lovelock terms correspond to the cases where all the 
vielbeins are the same. 

If the wedge products do not involve derivative terms, they 
are nonlinear potential terms for spin-2 helds 

•^pot — E^^^ A ■ ■ ■ A E^'^\ (3) 

which include the cosmological constant term and other in¬ 
teracting potentials for spin-2 helds E&lil in the vielbein 
formulation dB These potential terms are usually free of 
the Boulware-Deser (BD) ghost B For simplicity, they are 
denoted by de Rham-Gabadadze-Tolley (dRGT) terms B- 

The search for new derivative interactions was initiated in 
Cl where a new BD-ghost-free term in 4d was discovered 

(4) 

This cubic term can be thought of as a generalization of the 
perturbative Lovelock and dRGT terms and there are more 
possible terms in higher dimensions lfH . It was also conjec¬ 
tured that they should admit nonlinear completions in parallel 
to the dRGT terms in the context of massive gravity Cl- 

Some of the novel kinetic terms in ([1]) are nonlinear, multi¬ 
gravity completions of this cubic term in terms of differential 


^ In 1970, J. Wess and B. Zumino already proposed to use vielbeins, rather 
than metrics, as the building blocks of the low energy effective potentials 
for spin-2 fields 01. However, it is still not fully clear why the vielbein 
formulation is more efficient in eliminating the Boulware-Deser ghost. 

In a concrete model, the BD ghost may be present even if the building 
blocks themselves are free of it. For example, the BD ghost is excited if 
loop-type interactions are introduced in the metric formulation Ea. Anal¬ 
ogously, we find constraints on novel derivative terms. 

^ We refer to fT3l for reviews of this subject. 
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fonns@. If we consider two spin-2 fields, impose the symmet¬ 
ric condition iS and fix the second metric to Minkowski, 
the novel kinetic terms reduce to the two-derivative terms 
proposed in ifl^ . which are nonlinear derivative terms for 
a massive graviton around a Minkowski background. They 
can also be obtained by dimensionally deconstructing the 
Gauss-Bonnet term S. 

In the discussions of in and m . only one metric is 
dynamical and the Boulware-Deser ghost was shown to be 
present. In this work, we will not make the single dynami¬ 
cal metric assumption, and the conclusion is different con- 
cerning the fate of the Boulware-Deser ghost. 

In new kinetic interactions in 4d with a dynamical 
metric and a fixed Minkowski metric were investigated in de¬ 
tail. In the minisuperspace approximation, problematic N~‘^ 
terms were found in the Lagrangians and the corresponding 
Hamiltonians are nonlinear in the lapse function N. This 
indicates the secondary constraint from the time derivative of 
TTjv = 0 is an equation for N Q. Then the dangerous, sixth 
degree of freedom remains dynamical. In 4d massive gravity, 
the sixth degree of freedom is ghost-like. It plagues a generic 
nonlinear completion of Fierz-Pauli theory and is known as 
the Boulware-Deser (BD) ghost 0. This ghost-like degree 
of freedom is eliminated in nonlinear massive gravity @] 
thanks to the existence of a secondary constraint and an as¬ 
sociated tertiary constraint 1^ . Besides the minisuperspace 
discussion, an impressive no-go theorem was established in 
Cl on new kinetic interactions for single dynamical met¬ 
ric models that are Lorentz-invariant and free of the BD ghost. 

Inspired by the successful extensions of massive gravity 
H to bi-gravity H and multi-gravity d, we would like to 
examine the bi-gravity models involving the novel derivative 
terms O- Given that pathologies were found in single dynam¬ 
ical metric models, it is very likely that the bi-gravity theories 
are sick as well. In fact, more recently, bi-gravity models with 
new kinetic interactions were studied in the first order for¬ 
mulation d, where negative results were presented again0. 
Other obstructions were discussed in 1^ as well. 

Contrary to the single dynamical metric models, our 
analysis in section |III] show that the sixth degree of freedom 
could be absent in some bi-gravity models constructed from 
©. But the price to pay is that one of the linearized kinetic 
terms has a wrong sign or at least one of them vanishes. In 
the former case, we encounter spin-2 ghosts, which can lead 
to tree-level non-unitarity upon quantization. In the latter 


^ When the first version of this work appeared in arxiv, we were informed 
that similar nonlinear completions were studied before by K. Hinterbichler 
and R. A. Rosen. 

^ The absence of an additional constraint was already found in 0, which 
should be present after a change of variables. 

* By an abuse of terminology, the sixth degree of freedom of a massless 
spin-2 field in 4d is also denoted by the Boulware-Deser (BD) ghost. 

^ To avoid confusion, let us emphasize that we consider second order for¬ 
mulation in this work, so spin connections are not independent. In other 
words, we assume the torsion-free condition. 


case, the bi-gravity theories are strongly coupled due to the 
absence of kinetic terms. 

The presence of a spin-2 ghost is a well-known feature of 
a generic model of higher derivative gravity. In this work, 
we propose a new viewpoint that, when the couplings to mat¬ 
ter are not introduced, a large type of bi-gravity models are 
equivalent to higher derivative gravity without Ostrogradsky’s 
scalar ghost. They include Weyl gravity, 3d New Massive 
Gravity and some of their generalizations. 

At the classical level, a spin-0 ghost is more dangerous 
than a spin-2 ghost. Usually, the Hamiltonian of a scalar 
ghost is unbounded from below, while that of a massless spin- 
2 field may simply vanish. In this sense, it is more cru¬ 
cial to eliminate the Boulware-Deser ghost. Furthermore, the 
Boulware-Deser ghost should be removed if a massive spin-2 
field contains a correct number of dynamical degrees of free¬ 
dom, which is at most 5 in four dimensions. 

Upon quantization, a spin-2 ghost will lead to tree-level 
non-unitarity when coupled to matter. Let us remind that a 
notorious problem in quantizing gravity is that the Einstein- 
Hilbert action is non-renormalizable HH], which is very 
different from the other fundamental forces. By introducing 
higher order curvature terms (thus unitarity is sacrificed), one 
can improve the short-distance behaviour of the propagators 
and obtain a perturbative renormalizable theory for gravity 
0. Roughly speaking, the improved high energy behaviour 
are due to the relatively negative contributions from the ghost 
modes, which is analogous to the role of heavy ghost-like 
modes in Pauli-Villars regularization and superpartners in 
supersymmetric theories. 

From a different perspective, we can consider a metric as 
an effective description of some microscopic physics. In an 
effective field theory (EFT) of gravity, higher order curvature 
terms are expected in a low energy expansion, because they 
are compatible with the symmetries Even if gravity it¬ 
self is not quantized, they can be generated by quantum cor¬ 
rections from the matter. Therefore, some bi-gravity models 
with novel kinetic terms belong to a special class of effective 
field theories of gravity where Ostrogradsky’s scalar ghost is 
removed. 

The cut-off scale of an EET of gravity is usually associated 
with the Planck mass. However, if some of the high order 
curvature terms have unnaturally large coefficients, then the 
cut-off scale will be lowered. By eliminating the BD ghost, 
we could increase the cut-off scale set by large high order 
curvature corrections. 

This paper is organized as follows. In section [III we give 
the precise formulation of the novel kinetic terms for 4d bi¬ 
gravity models. In section [HIl we carry out a minisuperspace 
analysis to identify the candidate theories that are free of the 
dangerous, sixth degrees of freedom. In section we per¬ 
form a field redefinition to obtain more explicit expressions of 
the novel kinetic terms. In section |V] we focus on the novel 
kinetic term £ 2 '“" and discuss its relation to higher curvature 
gravity. In section IVTl we linearize the nonlinear models and 
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diagonalize the quadratic actions. In section IVlIl the issue of 
spin-2 ghost is discussed. In section rVlIIl we examine the 
cut-off scales of the effective field theories of gravity involv¬ 
ing novel kinetic terms. In sectionH^ we explain how to ob¬ 
tain the symmetric condition from the equations of motion. In 
section |X] higher derivative generalizations are discussed and 
a general argument for the absence of the BD ghost is pre¬ 
sented. In section IXIl we summarize our results and discuss 
their implications. 

II. NOVEL KINETIC TERMS EOR BI-GRAVITY 

To be more specific, we mainly consider four dimensional 
spacetime and models with two vielbeins/metrics. There are 


six possible nonlinear kinetic terms 

=R{E) ^E ^E, (5) 

£ 2 '^“ = R{E) AEAF, (6) 

£ 3 '^“ = RiE) AEAF, (7) 

£ 4 ^“ = R(F) AEAF, ( 8 ) 

£ 5 !^” = R{F) AFAE, (9) 

=R{F) AE AE, (10) 

where R{E), R{F) are curvature two-forms 

R{E) = du}^ +u}^ Auj^, (11) 

R{F) =duj^ +uj^ Auj^. (12) 

Both E and F are dynamical vielbeins. and u!^ are the 
spin connections compatible with E and F respectively 

D^E = dE + u;^ AE = 0, (13) 

D^F = dF + uj^ AF = 0. (14) 


Notice that one of the vielbeins in £ 2 *^“ and £ 4 '^"’ are La¬ 
grange multipliers. Another interesting observation is that 
£ 3 ’“" and can be thought of as the Palatini formulation 
of the Einstein-Hilbert term, where the spin-connections are 
expressed in terms of the associated vielbeins 0 
We also have five potential terms 


£lPOt ^ 

= E A 

E A 

E 

A 

E, 

(15) 

£^pOt ^ 

= EA 

EA 

E 

A 

F, 

(16) 

£3Pot ^ 

= EA 

EA 

F 

A 

F, 

(17) 

£4^°' = 

= EA 

E A 

F 

A 

F, 

(18) 

£5'’°' = 

= EA 

FA 

F 

A 

F. 

(19) 


When we discuss other dimensions, the subscripts n in £„ 
means that the number of F vielbeins in a wedge product is 


*** The difference is that varying the action with respect to the vielbein in 
the curvature two-form will give rise to second order equations. “Torsion- 
free condition" is not the only solution, so E and F are not necessarily 
proportional to each other. 


(n — 1). The two vielbeins E^^ and F^^ are related to two 
metrics and 

9^. = E^^E,,^gAB, Uu = F^^F,^r,AB. ( 20 ) 

The nonlinear kinetic terms £1 and £4 are the standard 
Einstein-Hilbert kinetic terms, while £ 2 , £ 3 , £ 5 , £5 are 
novel kinetic terms for spin -2 fields. 

To simplify our notation, the Levi-Civita symbol cabcd is 
not written explicitly in a wedge product. In terms of the com¬ 
ponents, the bi-gravity kinetic terms (ISlfTOli and the bi-gravity 
potential terms (I15II19I) are 


£ 1 “" = S^i^co R{E)^,^^Ep^E^°d^x, 

(21) 

£ 2 “" = R{E)p,,^^Ep^F^^d^x, 

(22) 

= STJ’co R{E)p,^^Fp^F„^d^x, 

(23) 

£ 4 '^” = 

(24) 

£ 5 “" = R{F)p,^^Fp^E^^d^x, 

(25) 

£ 5 “" R{F)p,^^Ep^E,^d^x, 

(26) 

£iP°‘ =S>17c% Ep^E,^Ep^E,^d'^x, 

(27) 

£ 2 ^°' = 6Zfco Ep^E.^Ep^F^^d^x, 

(28) 

£ 3 ”°' = SZ'’co Ep^E^^Fp^F^^d'^x, 

(29) 

= SZfJn Ep^F,,^Fp^F^°d^x, 

(30) 

£ 5 ^°' = S> 17 cd Fp^F,^Fp^F,°d^x, 

(31) 


A Fi 

where are the components of the curvature two- 

forms. The antisymmetric Kronecker delta or the generalized 
Kronecker delta is an antisymmetric product of the Kronecker 
deltas 

(32) 

where the antisymmetrization [... ] is not normalized. 
The Planck mass is not written explicitly. Greek letters 
fi, V, p, a,... indicate external spacetime indices and capital 
Latin letters A, B, C, D,... denote internal Lorentz indices. 

To minimalize the numbers of dynamical degrees of free¬ 
dom, we impose the symmetric condition or the Deser-van 
Nieuwenhuizen condition id. In section llXl we discuss how 
the symmetric condition originates in equations of motion. 


III. MINISUPERSPACE ANALYSIS 

In this section, we study the minisuperspace approximation 
of the bi-gravity models constructed from (ISlfTOl). The minisu¬ 
perspace analysis is a simple test of the ghost-free condition 
before investigating the fully nonlinear structure. Although it 
is not sufficient to prove healthiness, it is very efficient in de¬ 
tecting pathologies. Eor example, it was used in (idil to show 
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that loop-type interactions of multi-gravity in the metric for¬ 
mulation can excite the BD ghost. 

InS, new kinetic interactions with single dynamical met¬ 
ric were ruled out, because their Hamiltonians are not lin¬ 
ear in the lapse function in the minisuperspace approxima¬ 
tion, which indicates the presence of the BD-ghost. Despite 
the failure of single dynamical metric models, a large class of 
bi-gravity theories do satisfy the criterion that the minisuper¬ 
space Hamiltonians are linear in lapse functions, as we discuss 
below. 

Now we start the minisuperspace analysis. The two metrics 
in the minisuperspace approximation are diagonal 

dsi = dx^dx" = —{Ni^dt^ + e^^dx'^, (33) 

dsl = dx^dx’' = — (iV2)^df^ -f e^^dx^. (34) 

where the metric components are functions of time 


Pnx = Pn2 = 0; (42) 

where the second line are primary constraints. Then we derive 
the minisuperspace Hamiltonian by the Legendre transform 

P-mua = -^Pa + BPb — i^niini- (43) 

Using the relations between momenta and velocities, one 
can express the Hamiltonian T-Lmmi in terms of 

iV2, A, B, Pa, Pb- (44) 

The explicit expression of T-Lmim is a fraction 

Hmini = (45) 

4/xd 

where the numerator Pl^ and the denominator Pl^ are 


N^ = N^{t), N2=N2{t), A = A{t), 

The corresponding symmetric vielbeins are 


B = B{t). 

(35) 


E,. 


(Ni 0 \ 

'vO 


F,, 



Let us consider a linear combination of 
the Lagrangian reads 



(36) 

.., CfP'', so 


£ = ai £1“" + 02 £2"^“ + 03 £3'^" 

+6i£4“'’ + 62£5'‘“ + fc3£6“". (37) 


In the minisuperspace approximation, it becomes 

£mini = ai^{Afe^^ P bi^{Bfe^^ 
iVl 1^2 

+ 2ii? - g2A+B 

J-n ^ foAR ^2 12 A-b\ A+2B 

+a,— [2AB-—Ae Je 

where some normalization factors are inserted to simplify the 
expression of £niini- Time derivatives are denoted by dots 


Bn = ( 636 ^ + 626 ^) e- 2 ^ Pi (Afi)3 
-|-(a2e^ -I- 036-®) (^¥2)^ 

+2(636^ + 626^) PaPb { NifN 2 

+2(026^ + 036^) PaPb Ni { N 2 f 

-{b2e^ + bie^) Pi {NifN2 

-(aie^ + a2e^)e-2^P|7Vi(7V2)2, (46) 

B-d = (^ 36 ^ + ^ 26 ^) [(oi + b^e^ + (02 + ^ 2 ) 0 ^] (TVi)^ 
-|-(a2e'^ -I- 036-®) [(02 + 62)6^ + (03 -I- bi ) e ^] (iV2)^ 
- [(0162 - 0263 ) 6 ^"^ + ( 02^1 - 0362 ) 6 ^'® 

+(0161 - 0363 ) 6 ^+^]7Vi7V2. (47) 

We require "Hmini to be linear in Ni and N 2 , so Ni and N 2 
are Lagrange multipliers. Then the secondary constraints 


PN\ — ^ 

^0, 

(48) 

Pn2 — ^ 

s 0 

(49) 

are equations for the canonical variables and could remove the 
scalar modes related to the BD ghost[3- The Poisson bracket 
is defined as 

(“'«= L {s,ap- 

q=A,B,Ni,N 2 ^ ^ ® 

da 

dPq dq J 

(50) 


-A = A, 


-B = B. 


(39) 


dt dt 

Integration by parts is performed in order to eliminate 

Ni, N 2 , B. (40) 


The conjugate momenta can be derived from the minisuper¬ 
space Lagrangian 


The numerator Ha and the denominator are polynomials 
of degree three and two in Ni and 7 V 2 . To satisfy the require¬ 
ment that lapse functions are Lagrange multipliers, Hd should 
be a factor of Ha- This is true when only one of the three 
monomials in Hd has non-zero coefficient, which indicates 
two classes of bi-gravity models: 


Pa 


dCa 


dA 


Pb = 


dCa 


dB 


11 


(41) 


Strictly speaking, we should use the total Hamiltonian which contains the 
primary constraints to compute the time derivative. 
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• The first class is 

02 = as = 62 = 63 = 0, (51) 

where the Lagrangian contains two Einstein-Hilbert 
terms and the minisuperspace Hamiltonian reads 

= (52) 

One can introduce the cosmological constant terms 
(1151119b . The interactions between the two metrics are 
through the nonlinear potential terms dlblllTIfTST l. The 
minisuperspace Hamiltonian is still linear in the lapse 
functions. 

Since we have two Planck masses in front of two 
Einstein-Hilbert terms, one can take the limit where one 
of them goes to infinity. In this decoupling limit, a bi¬ 
gravity model reduces to that of single dynamical met¬ 
ric with a fixed metric. 


We can introduce the potential terms £ 1 ^°',..., £ 5 ^°*. 
The lapse functions will still be Lagrange multipliers in 
the Hamiltonians. 

Erom the holographic point of view, the diagonal dif- 
feomorphism invariance is fundamental in bi-gravity 
and multi-gravity theories 12^ . In the context of the 
AdS/CET correspondence, a conformal field theory 
with conserved stress tensor is dual to a diffeomorphism 
invariant theory. The massive gravitons (or more gen¬ 
eral the spin -2 fields without gauge invariance) corre¬ 
spond to spin-2 operators that are not conserved. Erom 
this perspective, a massive gravity theory should always 
admit enhancement to a diffeomorphism invariant the¬ 
ory by turning on the conserved stress tensor in the 
boundary field theory. However, the converse is less 
justified. A bi-gravity or a multi-gravity theory may 
not have a decoupling limit that breaks the diagonal dif¬ 
feomorphism symmetry, which indicates the conserved 
stress tensor decouples. 


• The second class is 


ai = 02 = 03 = 0 

(53) 

b\ = b2 = bz = 0 . 

(54) 


The bi-gravity models in the second class contain at 
most one Einstein-Hilbert kinetic term (no EH term if 
ai = bi = 0 ), so one can not take the decoupling limit 
that fixes one of the metrics [^[EIH Therefore, they are 
not ruled out by the no-go theorems in ifT^ . We focus 
on the case of ( |54] | in the discussions below. 

The minisuperspace Hamiltonian with (l54l i reads 


Km = 


+N 2 


2(02 -b 03 eB) 

(oi -b 02 
4(02 -b 03 eBy 

2 ' 

4(026^ 4-036-^)^^ 



(56) 


where we assume 02 and 03 are not zero at the same 
time. 


There exists another single metric limit 

E — F = H/X, X —>■ 00, ( 55 ) 

where the novel kinetic terms reduce to the Einstein-Hilbert term. 

The failure of obtaining nonlinear partially massless gravity from a consis¬ 
tent truncation of conformal gravity (H is related to the absence of this 
decoupling limit. 

The fact that only certain combinations of kinetic terms are allowed is 
analogous to the scalar-tensor theories discussed f^ . where the degen¬ 
eracy conditions can break down for some combinations of degenerate La- 
grangians. 


According to the minisuperspace Hamiltonians, there are 
two classes of bi-gravity models that are potentially free of 
the Boulware-Deser ghost. The first class of bi-gravity the¬ 
ories were proposed in by promoting the fixed metric in 
consistent non-linear massive gravity |0] . Only the standard 
Einstein-Hilbert kinetic terms are used. The absence of the 
BD ghost was proved in 

The second class of bi-gravity models (l54l) are in a different 
region of the parameter space where novel kinetic terms are 
used. Let us investigate them in more detail. We can compute 
the Poisson bracket of the two constraints associated with Ni 
and N 2 


oKm dm 


dN, 


dN2 


„ „- 2 B p 2 

03 e fg 


oi ^Pb 


8(02 -b 03 eBy 

—03 Pa — 02 { Pa — Pb ) 


(57) 


If 02 ^ 0 and 03 = 0, the Poisson bracket (l57l i vanishes. It 
seems that the two constraints are first class constraints, which 
could be associated with two sets of gauge symmetries. If 
02 0 and 03 ^ 0 , an independent constraint is generated 

by the stability of secondary constraints, which involves cu¬ 
bic momentum terms according to ( fSTT i. This constraint elim¬ 
inates one more dynamical variable, which signals the com¬ 
plete absence of the sixth degree of freedom. 

Note that it is not justified to take the minisuperspace ap¬ 
proximation before computing Poisson brackets. But this 
simple computation does capture the essential features of the 
Hamiltonian structure: 


1. when 03 = 0, the Poisson brackets of secondary con¬ 
straints vanish on the constraint surface; 

2 . when 03 7 ^ 0 , an independent tertiary constraint is gen¬ 
erated, rendering absent the sixth degree of freedom. 
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At this point, it is not clear whether the minisuperspace re¬ 
sults can be extended to the full theories. To verify this, we 
need to examine the Hamiltonian structure of the full theories, 
which is highly technical. 

In the first class of bi-gravity models (fSTl i. the kinetic terms 
are the standard Einstein-Hilbert terms. The secondary con¬ 
straints are those in general relativity supplemented by the 
contribution from the potential terms (llblllTIfTSl l. which do 
not involve momenta and spatial derivatives. To compute the 
Poisson brackets of the constraints, one can use Dirac’s hy¬ 
persurface deformation algebra. 

However, in the second class of bi-gravity models (l54l i. the 
kinetic part of the Lagrangians are modified by the novel ki¬ 
netic terms. To obtain the Hamiltonians already requires some 
work. The constraints have more involved dependence on mo¬ 
menta and spatial derivative terms. Dirac’s algebra is not ap¬ 
plicable. The computations of constraint brackets are consid¬ 
erably more challenging. We leave the technical analysis of 
the Hamiltonian structure to a separate work ii in which we 
verify that the sixth degrees of freedom are indeed eliminated, 
and the case 02 0 , 03 = 0 does describe two interacting 

massless spin -2 fields. 


IV. FIELD REDEFINITIONS 

Let us derive the explicit expressions of the novel kinetic 
terms. It is difficult to achieve this step directly because the 
components of a curvature two-form are complicated func¬ 
tions of the associated vielbein. To circumvent this difficulty, 
we make use of a mathematical identity for tensors in d di¬ 
mensions!^ 

det(£;)T';i-‘^‘'(E;-i);.J''' ... (58) 

where the antisymmetrized product of E~^ gives det(£'“^) 
and cancel det(£') out. For example, in two dimensions we 
have 

rpOl _ ^10 

= det{E)T°^[{E-^)o°iE-^)^^ - iE-^)o\E-^)i°] 

+ det{E)T^°[iE-^),°{E-^)o^ - {E-^),\e-^)o°]. 

We notice the minisuperspace Lagrangian (l3^ contains the 
ratio of two lapse functions. From the mathematical identity 
(l58l l. it is natural to introduce a new tensor field as the “ratio” 
of two vielbeins 

= (59) 

where p and ly are external spacetime indices and A is an in¬ 
ternal Lorentz index. 

Using the identity (l58l l. the novel kinetic terms become 
£ 3 “" = d*x . (60) 


By an abuse of notation, the local Lorentz indices are denoted by Greek 
letters as well. 


After simple manipulations, we have 

£ 2 *"” = (-)^ (61) 

R{9)iJ.u^^''epPea'^^ d'^x, (62) 

where the antisymmetrization [... ] is not normalized. The 
normalization factors in dMT l and ( l62l i are made precise. These 
choices simplify the expressions below. More explicitly, £ 2 ’'“ 
and £ 3 ’“" are functions of the metric g^i, and the new tensor 
field ep,’' 


£ 2 '^" = 


(63) 


>^ 3 '^” = R^u'^'^ep'^e,'' - 2i?/( 




Pp pt 




d^x. 


(64) 


In terms of the new tensor field e^v, the potential terms (flSl 
[Toll are 


= (65) 

£2P°‘ = v/^ep^ (66) 

Cr^ = ^gep^^^e,'^\ (67) 

Cr' = ^f^ep^^^e,'^epP\ ( 68 ) 

£ 5 ?°' = y^ppl'^ep'^ep^e^'"!. (69) 


We can lower and raise the indices of by the metric g^i, 
and its inverse g^^'^. For example, we have 

epv = epPgpv = Fp^E^^gAB- (70) 


The symmetric condition then translates into 

ep^ = Fp^E^^gAB = F^^Ep^gAB = e„p (71) 


and we have 


e^p eP„ = g^°'eai3g^^ep„ = fa„. (72) 

When Bpi, is symmetric, e^i, is the square root of g~^f with 
fiJ.i' defined in I®. 


V. £2*“" & QUADRATIC CURVATURE GRAVITY 

In this section, we focus on the novel kinetic term £ 2 '“" and 
explain its connections to the models of quadratic curvature 
gravity. 


A. Gauge symmetries of £2'^'" 

As anticipated in the minisuperspace analysis, the case of 
03 = 0 , 02 0 seems to be related to bi-gravity models 

with two sets of gauge symmetries. Let us emphasis that the 
coefficient oi of the Einstein-Hilbert term is not fixed. 

From the explicit form of £ 2 '“", we can see linear combi¬ 
nation of £ 1 '^” and £ 2 ^^” are invariant under several gauge 
transformations: 
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• Standard diffeomorphism invariance 


B. New Massive Gravity 


(73) 

where is a four-vector. This symmetry is expected in 
a covariant bi-gravity theory. 

• Additional “diffeomorphism invariance" 

(74) 

where is a four-vector. The Lagrangians are invariant 
because £ 2 *^“ is the product of Einstein tensor and the 
new tensor After integrating by parts, the covariant 
derivative acts on the Einstein tensor and the change in 
the Lagrangians vanish due to the second Bianchi iden¬ 
tity. 

• Local Lorentz invariance 


One may suspect the bi-gravity models with 03 = 0 can 
be transformed into two Einstein-Hilbert terms after some 
field redefinitions. It is not clear what field redefinition can 
make this connection. Nevertheless, it was shown in ll^ that 
£ 2 '“" in (l 6 ^ can be obtained by taking a scaling limit of two 
Einstein-Hilbert terms 


V^7?(/) - yf^R{g) 




(80) 


where is defined as 




( 81 ) 


(75) with 


Eor an infinitesimal transformation, we have 

5E^^ = UJB^ = iOB^ F^^, (76) 

where 

Ab^ = Sb^+ujb^+0{uj^). (77) 

Erom the definition (|59] |, we know is invariant un¬ 
der a diagonal local Lorentz transformation, so the La¬ 
grangians are invariant as well. 

• Additional “local Lorentz invariance" 

- ^I^fl (78) 

or in the infinitesimal form 

JF / = uj'b^ F /, u;'^b = -^'ba- ( 79 ) 

The antisymmetric part of is projected out by the 
symmetric Einstein tensor, so the Lagrangians are in¬ 
variant under a change in the antisymmetric part of e^i/. 

Since the antisymmetric part drops out when 03 = 0, 
we can identify the symmetric part of e^i, with the 
square root of a metric product E 

These gauge symmetries persist even if we turn on the po¬ 
tential terms £iP°‘ and £ 2 P°^ which are related to the cosmo¬ 
logical constant. However, the “additional" gauge symmetries 
will be broken when £3'"", £3^°', £4^°^, £5^°^ are introduced. 

It should be noted that the additional “invariances" are not 
the precise transformations according to their names, but they 
do have the same amount and similar forms of gauge symme¬ 
tries as indicated. To be more precise, if we substitute Se^i, on 
the left hand side of IT^ with Sg^i, and SFp^ on the left hand 
side of ( l79l l with SE^^, they become the off-diagonal gauge 
transformations. 


A —>■ 00 . (82) 

In fact, £ 2 '^“ in the form of (l 6 ^ already appeared in the 
auxiliary field representation of 3d New Massive Gravity ll^ 
in, which is a theory of quadratic curvature gravity. 

In our language, the Lagrangian of 3d New Massive Gravity 
reads 0 

EnMG = Cr£EH + £ 2 *^” + Ci£iP°‘ + C3£3P°‘, (83) 


£eh = R[E) A F, (84) 

£ 2 “" = R{E) A F, (85) 

£iP°‘ = EAEAE, ( 86 ) 

£ 3 P°‘ = F a F a F, (87) 


where cr = ±1 is the sign of the Einstein-Hilbert term, ci is 
proportional to the cosmological constant, C 3 corresponds to 
the mass squared m^. The usual auxiliary field is identified 
with the second tensor field 

= epPgpv = Fp^E^^rjAB- ( 88 ) 

Note that the symmetric condition 

(89) 

is imposed dynamically by the equations of motion 0 - 

Therefore, New Massive Gravity is an example of 3d bi¬ 
gravity models in the second class. It is known that 3d New 
Massive Gravity do not contain the Boulware-Deser ghost 
lf33l . which furnishes evidence that the second class of bi¬ 
gravity models are free of the BD ghost. 


Note that the metric is different from defined in |20j and they 
coincide only when the antisymmetric part of e^u vanishes. 


In this representation, we can see one of the BD-ghost-free potentials al¬ 
ready appeared in the context of New Massive Gravity. 

In section E] we discuss how to generalize this example of dynamical 
symmetric condition. 
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A straightforward generalization of 3d New Massive Grav¬ 
ity is to introduce other potential terms £ 2 '’°': £4 

£2P°‘ = EAEAF, (90) 

£4P°‘ = FAFAF. (91) 

Note that, if £ 4 P°'^ is considered, the Lagrangian does not re¬ 
duce to that of quadratic curvature gravity when is inte¬ 
grated out. Instead, it contains infinitely many higher order 
curvature corrections. 

In 3d New Massive Gravity, the Einstein-Hilbert term is al¬ 
ways present. One of the reasons may be that the second spin- 
2 field is usually considered to be an auxiliary field, which 
seems to have no dynamics. We want to emphasize that £ 2 '“" 
is a kinetic term as well, and the use of the Einstein-Hilbert 
term is not necessary 0- Therefore, the Einstein-Hilbert term 
could be absent, then the 3d bi-gravity Lagrangian reads 

4 

£ = 02 £ 2 *^" + ^ C, £,P°', 02 ^ 0. (92) 

i=l 

The number of dynamical degrees of freedom should be the 
same as that of New Massive Gravity and the Boulware-Deser 
ghost should not be not propagating. The cases without the 
Einstein-Hilbert term are related to the generalized NMG in 
by a field redefinition 

^ T cp/iiv. (93) 

where c depends on the coefficients of £eh and £ 2 ^^”. 

In 3d, £ 2 '^“ is the only novel kinetic term from ([T]i due to 
the limited number of spacetime indices. Since a massless 
graviton in 3d has no dynamical degree of freedom, we can 
choose 02 such that the kinetic term of the massive graviton 
has a correct sign. Then ( |9^ is a unitary theory of 3d massive 
gravity. The special case of ci = C 2 = C 4 = 0 was discussed 
in and that of ci 7 ^ 0 , C 2 = C 4 = 0 in j^ . 

C. Critical gravity 

There exists a continuous family of critical points in 
the parameter space of 3d New Massive Gravity dS^ and its 
higher dimensional generalization , 

C = R{E)AEA---AiE + E) 

+E A-■■ A E A {AE A E+ m^ E A E), (94) 

which are known as critical gravity. At these critical points, 
the cosmological constant A is proportional to the mass 


To eliminate the second order time-derivative terms due to the curvature 
tensor in ^ 2 *“", we need to supplement the action by boundary terms anal¬ 
ogous to the York-Gibbons-Hawking term. Then € 2 ^'" generates a time 
derivative term dtefiv in the Lagrangian, so both g/^, and have dy¬ 
namical degrees of freedom. Furthermore, if we expand the Lagrangian 
around a Minkowski background and diagonalize the quadratic kinetic 
terms, £ 2 ^“ will give rise to two linearized Einstein-Hilbert terms. 


squared m? with dimension-dependent coefficients. Integrat¬ 
ing out the auxiliary field, the linearized 4th order equation of 
motions contains two massless spin -2 modes PI 
Here we want to point out that 

£ 2 “" = R{E) A E (95) 

and its higher dimensional version 

£ 2 '^” = R{E) AEA---AEAE (96) 

can be thought of as a special limit of critical gravity in the 
bi-gravity formulation PI. Note that in this limit, can not 
be integrated out because it is a Lagrange multiplier. 

It is shown in the Hamiltonian analysis of 11^ that £ 2 *^” in 
4d has two sets of first-class constraints, corresponding to two 
sets of gauge symmetries. More general, the Lagrangian 

£ = ai £1'^” + 02 £2'^” + Cl £iP°‘ + C2 £2P°‘ ( 97 ) 

with 

02 ^ 0 (98) 

are interacting theories of two massless gravitons in various 
dimensions (d > 2 ), where Ci and C 2 are related to the cosmo¬ 
logical constantly. 

The same gauge symmetries are also realized in higher 
derivative counterparts of the two-derivative term £ 2 '“" 

R{E) A • • • A R{E) AEA---AEAE, (99) 

where one of the E-vielbeins in the Lovelock terms IH] is 
replaced by an E-vielbein. The additional symmetries are 
due to the fact that Lovelock tensors are both symmetric and 
divergence-free. 

Along the line of the second class bi-gravity theories (l54l) . 
we propose a general Lagrangian describing two interacting, 
massless, gauge invariant gravitons, which is a linear combi¬ 
nation of Lovelock terms, the novel derivative terms (l96l l99b 
and two potential terms 

= E A--- A E, £ 2 P°‘ = £A---A£AF, (100) 

where at least one of the novel derivative terms is presentPI. 


However, the total number of dynamical degrees of freedom should be the 
same as that of 1 massless and 1 massive gravitons. The second massless 
graviton seems to be an artifact of the linearized equation of motion at 
the critical points, as there is no symmetry enhancement. For example, 
logarithmic modes ai‘e allowed if we do not assume the Brown-Henneaux 
boundary conditions (3^ . They were claimed to be dual to Logarithmic 
Conformal Field Theories f^l4Qlj43l . 

The Einstein-Hilbert term is absorbed into by redefining F. 

There is a no-go theorem for interacting theories of massless, gauge invaii- 
ant, spin-2 fields if the Lagrangian has at most two derivatives f43ll . This 
is not in contradiction to the present work, because one of the linearized 
kinetic terms has a wrong sign, which violates one of the assumptions in 
the no-go theorem. The details of the linearized actions are discussed in 
section |vi] A recent construction of color-decorated gravity B evades 
this no-go theorem by including extra fields. 

Deforming these massless models by other potential terms with two viel- 
beins, one obtains the generalizations of New Massive Gravity proposed in 

ill]. 
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D. Weyl gravity 


Weyl gravity is a well-known theory of quadratic curvature 
gravity in 4d, which is both diffeomorphism and conformal 
invariant. Interesti^ly, the Lagrangian of Weyl gravity can 
be reformulated asFi 


£ 2 '^“ and £ 3 ^^” are the novel nonlinear kinetic terms defined 
in (EHISal). 

Let us expand the metric field and the symmetric tensor 
field around the Minkowski background 

(104) 


Cweyi = R{E) AE A F + E A E AF A F, (101) 
where F has dimension 2. 

Then the absence of Ostrogradsky’s scalar ghost is trans¬ 
lated into the absence of the BD ghost. Ostrogradsky’s spin-2 
ghost in the four-derivative formulation now becomes a basic 
spin -2 ghost due to a wrong sign kinetic term. 

In this representation, Weyl gravity is built from a novel ki¬ 
netic term and a dRGT potential term. Despite the presence of 
a spin-2 ghost, Weyl gravity is the first example of nonlinear 
completions of Fierz-Pauli massive gravity that are free of the 
BD ghost, where the nonlinear theory was proposed 20 years 
before the linear one. 

Furthermore, Weyl gravity is a special bi-gravity model 
in the second class ( |54| ) with an emergent gauge symmetry 
(conformal symmetry). This gauge symmetry is a nonlinear 
completion of the additional gauge symmetry of a massive 
spin-2 field around de-Sitter background at the partially mass¬ 
less point iz^ . To make this connection more clear, we lin¬ 
earize ( IIOII ) around de-Sitter background and diagonalize the 
quadratic Lagrangian in the next section. 


VI. LINEARIZED LAGRANGIANS 

In the previous section, we show that some of the bi-gravity 
models with novel kinetic terms are equivalent to higher 
derivative gravity models. It is well known that higher deriva¬ 
tive gravity models usually contain spin -2 ghosts, which could 
lead to the problem of non-unitarity. In this section, we derive 
the quadratic actions of the novel kinetic terms and examine 
whether this is a general feature of the bi-gravity theories in 
the second class (l54l) . 


A. Minkowski background 

Consider a bi-gravity model in 4d whose Lagrangian reads 
£ = oi £eh + 02 £ 2 '^™ + 03 (102) 

where £eh is the Einstein-Hilbert kinetic term 

Eeh = y/^R{g), (103) 


By redefining F = F' + E, one have the auxiliary field reformulation 
of Weyl gravity with an Einstein-Hilbert term and a cosmological constant 
term, where the second spin-2 field can be thought of as a matter field 
couples to the Einstein tensor. 


where we assume there is no numerical factors in front of rj^^. 
These factors can always be absorbed into ai, 02 , as by re- 
defining and 

Note that simply vanishes in a different kind of back¬ 
ground solutions. The two kinds of background solutions are 
related by a redefinition of F^^^ 

F' = F + E. (106) 


To the quadratic order, the linearized Lagrangian reads 
£ = Cl + C 2 (107) 


where the first term is the linearized Einstein-Hilbert term and 
the coefficients are 

1 X 1 

Cl =-(-ai-I-02 - as), C2 = -a2-a3. (108) 

Now we can diagonalize the quadratic Lagrangian 

£ = Cl (h^^^d,d’'h/ - , (109) 


where we assume ci 0 and the diagonalized spin -2 fields 
are 


C 2 C 2 

“b 7i (110) 

Zci Zci 

If Cl = 0 and C 2 yb 0, then the first term in ( 1107b vanishes 
and the diagonalized Lagrangian is 


£ = 


C 2 , ( 111 ) 


and the diagonalized fields are 

hfii/ — Fpijj — — {Sg^ii! ). ( 112 ) 


The linearized Lagrangian after diagonalization is a linear 
combination of two linearized Einstein-Hilbert terms 0 ex¬ 
cept in some special cases. The diagonalized kinetic terms al¬ 
ways have opposite signs due to the absence of quadratic term 
of e^i^, which can be traced back to the absence of R{F). 


The diagonalized form is invariant under a continuous family of field re¬ 
definitions 


fi/rv = cosh{9)h'^^ + smh(9)H'^^, 

= sinh(61)h),„ + cosh(d)H^^. (113) 
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The special cases are 

C2=0, or ^02=03, (114) 

then H^i, = 0 and the second diagonalized kinetic term van¬ 
ishes. The first kinetic term has a right sign when 

Cl = i(-ai + 03 ) < 0. (115) 

A more extreme case is 

Cl = C 2 = 0, or 01 = ^02=03, (116) 


We assume the background value of vanishes, which 
is not necessary. However, if the background solution of 
is proportional to the background metric then we can al¬ 
ways set it to zero by a shift in 

To avoid strong coupling problem, we require 

02 ^ 0 , ( 121 ) 

otherwise we should consider a different background solution. 

Around a background solution, the linear terms in the per¬ 
turbative Lagrangian should vanish, which indicates 

bi = —2ai A, 62 — tt 2 A. (122) 


then the linearized Lagrangian is empty. In these special 
cases, the bi-gravity models are strongly coupled due to the 
absence of some linearized kinetic terms. 

Note that the Lagrangians of these special cases can be 
schematically written as 

£ = (ai — as) RAE AE + as RA{F — E) A{F — E). (117) 

If we consider the other kind of background solutions where 
vanishes, both of the two linearized kinetic terms are 
present and there is no issue of strong coupling. So the prob¬ 
lem of strong coupling depends on the choice of background 
solutions. Ej The existence of these strongly coupled back¬ 
grounds is related to the presence of spin -2 ghosts. 

The potential terms ( I65H69I I can generate linear terms 
around a Minkowski background, which signals a wrong 
choice of vacuum. In the next subsection, we discuss the 
linearized actions around general maximally symmetric back¬ 
grounds. 


B. Constant curvature background 


Before the shift in F, the two equations correspond to the so¬ 
lution of the background metric g^i, and the ratio between two 
background spin -2 fields. 

The linearized Lagrangian is 


+V^C2 +,5e/[Vp, 

2aiA - aiA(^^Sg^^Sgi,'' - ^Sgf,''Sg^'"^ 

+63 {Se^^^Se,’' - (123) 

where the coefficients 


Cl = C 2 = -02 


(124) 


are simplified due to a shift in F. Total derivative terms are 
neglected. 

If oi 7 ^ 0, the diagonalized fields are 


Let us introduce nonlinear potential terms to the 4d La¬ 
grangian 


— bgfiu 

Ol 


TT _ 
fiv - 

ai 


(125) 


£ — Ol £eh + 02 £2'^” + 03 £3*^“ 

+ 6i£iP°‘ + 62£2P°‘ + ()3 
+ 64£4P°‘ + 65£5'’°‘, (118) 

where the potential terms are defined in (I65l l-(l69]l. 

The spin-2 fields and e^i, are expanded around a cos¬ 
mological background p^jy 

C/yiy — (119) 

where the background spacetime has constant curvature 

- A, _ _ _ ^ 

R/iiypa — ~^\9p.p9vcr 9pcr 9vp)- ( 120 ) 


the linearized Lagrangian becomes 
1 


£ = --oiv/^ V^]KP^ 

+ \aiy/^ +7T/[Vp, 


+aiA(- 


1 


\ 02 Z 


(126) 


where the constant term is neglected. The first four lines are 
the linearized Einstein-Hilbert terms with cosmological con¬ 
stant A around the background solutions. The last line is the 


If we consider backgrounds that vanish, then will not contribute 
to the linearized Lagrangian. We will encounter the strong coupling prob¬ 
lem when a 2 = 0. 


The definitions of fluctuating fields are modified accordingly. 
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Fierz-Pauli mass term for The coefficient of the mass¬ 
less spin-2 field is ai, while that of is — ai. One 
of them is a spin-2 ghost. The mass squared of H^i, is deter¬ 
mined by & 3 . 

If oi = 0, the diagonalized fields are 

T H = 2 (127) 


In addition, 63 should vanish in order to be consistent with 
our choice of background solution = 0 , so the mass terms 
vanish. The linearized Lagrangian is 


C = \a2V^ VP]KP) 


-\a2yfl +iT^"[Vp, V^]i7/) 


+\/—g 2 a 2 ^(^-hfj,P'h^'^ — 

-2a2A(ii7/i7," - 


, (128) 


where we keep the non-zero background value of e^j^. 

The fluctuations around the de-Sitter vacuum are 

( 134 ) 

Then we can expand the full action to the quadratic order 
-Cweyi = ^ “ Se^^deP'"') 


+ ^{SeSg - iSe^Jgn - ^{Sg^ - 4.6g^Jgn 
+lv^ {se^^^y,v''Sg/ + Se^''[yp, yp]Sg^p) 
~Ya'^ + (5fl/[Vp, ^P\5g,p) , 


(135) 


where total derivative terms are neglected, the covariant 
derivative V is compatible with the background metric g^^ 
and 


which corresponds to two interacting massless gravitons. 


^9 = Se = Scpp. (136) 


C. Linearized Weyl gravity 


In this subsection, we would like to discuss the linearized 
action of Weyl gravity ( llOlb around the de-Sitter background. 
As shown in , the conformal transformation in Weyl grav¬ 
ity can be recast into nonlinear partially massless transforma¬ 
tion for spin-2 matter field after some field redefinitions. So 
we expect after diagonalization the massive spin -2 field has 
linear partially massless gauge symmetry. 

The explicit expression of (IIOII) is 


■l^Weyl — 


kin 


.£3Pot 




-fiiy 


_Lp _ p 

where the equations of motion for is 


(129) 


1 


1 




(130) 


The nonlinear gauge symmetry transformations are 


The diagonalized fields are 

hpi/ = Hpiy = Sg^i/ (137) 

The quadratic Lagrangian in terms of H^i, reads 
Liweyl 

= [hp^PV,V''hpP'' + /i/[Vp, 

“4^ + iT/[Vp, V'^]i7,^) 

, (138) 

where we neglect the constant term. We can see ft. is a mass¬ 
less spin-2 field with a negative Planck mass, while H is a 
massive spin-2 field with a positive Planck mass. The signs of 
the kinetic terms are opposite. 

We can further examine the gauge symmetries at the lin¬ 
earized level: 


• conformal invariance 


• linearized conformal symmetry 


gfiiy ^ (1 “f 20)p^i/, Cpu y Gpu (131) 

• diffeomorphism invariance 

gfiu ^ 9tJ^i^ ^ (-flu + £^ep,. (132) 

Let us consider the de-Sitter background solution 

Hfiu = g^uy ~ 'q9pv^ (133) 


^gfiv —i" Sg^i, -I- 2(j)gp,y, 


(139) 

^(flU ^ ^(flU ^ fl^u4^^ 


(140) 

and 



6 - 





(141) 

f- ^ A_ \ 

/6 \ 


Hfiv Hfiu + ( -f '^9fiv j 

■ 

(142) 
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• linearized diffeomorphism symmetry 

+ V, (143) 

Se^i, —> Sefj,iy + — , (144) 

and 


sign of its kinetic term. However, there is only one new kinetic 
term in 3d, and these bi-gravity models were well investigated 
along the line of 3d New Massive Gravity. In particular, the 
3d version of Weyl gravity proposed in ll^ is an example 
of the conflict between diffeomorphism and conformal invari¬ 
ances in three dimensions. Partially massless symmetry has 
no nonlinear extension in 3d up to date. 


-f (145) 

Hf,. ^ (146) 


VII. SPIN-2 GHOST 


Therefore, the linear partially massless gauge transforma¬ 
tion 


H, 


fiiy 


H, 


fiiy 




(147) 


is a combination of linearized conformal and diffeomorphism 
transformations. 

Interestingly, only the massless spin-2 field h transforms 
under a change of coordinate. The massive spin-2 held can 
not transform because the Lagrangian of massive mode is not 
invariant. 

In this way, we provide a different perspective of Weyl 
gravity by using a novel kinetic term and a dRGT term. In 
this representation, one may understand better why unitary 
partially massless gravity in 4d is not found Pi . Along 
the line of dRGT massive gravity, there have been many recent 
investigations on nonlinear partially massless gravity 
In 4d, a promising candidate was identified in dRGT massive 
gravity ll4^ , which makes use of precisely the same potential 
term 


EAEAFAF, (148) 

but the kinetic term is assumed to be the Einstein-Hilbert term 
and the F vielbein is fixed to be de-Sitter. Partially mass¬ 
less gauge symmetry is only an artifact of the perturbative 
Lagrangian at low orders. We also confirm the suspicion in 
||4^ that a new kinetic term is required in order to extend the 
partially massless symmetry to the nonlinear level, which be¬ 
comes trivial from our bi-gravity reformulation of Weyl grav¬ 
ity. 

The fact that the novel kinetic terms have no non-trivial sin¬ 
gle dynamical metric limit indicates that we can not truncate 
Weyl gravity in a consistent manner to obtain a nonlinear the¬ 
ory of partially massless gravity with single dynamical metric 
and a fixed de-Sitter metric. 

In addition, partially massless gauge symmetry in 4d can 
be thought of as an emergent gauge symmetry of the 4d 
bi-gravity models constructed from novel kinetic terms and 
dRGT potential terms. It is tempting to consider the case of 
three dimensions, where a massless spin-2 held has no dynam¬ 
ical degree of freedom and we do not need to worry about the 


To be more precise, we require the off-shell action to be gauge invariant. 


In the previous section, we show that, generically, the two 
linearized kinetic terms have opposite signs. A kinetic term 
with a wrong sign usually results in an unbounded Hamil¬ 
tonian and non-unitarity. If we associate the Hamiltonian 
with the energy, then we could extract infinite energy from a 
system whose Hamiltonian is unbounded from below, which 
leads to classical instability. Non-unitarity in a quantum 
theory means unphysical negative probability. Therefore, 
when a kinetic term has a wrong sign, the corresponding 
degree of freedom is considered to be an unwanted “ghost". 

From the classical point of view, it is not clear which is the 
correct sign for a massless spin-2 kinetic term, as the Hamil¬ 
tonian simply vanishes on-shell. Naively, the Einstein-Hilbert 
term seems to have a wrong sign in the minisuperspace ap¬ 
proximation (El 

V^R{g)^ N>0, (149) 

but the Hamiltonian is still bounded because it vanishes. If we 
modify the sign of the Einstein-Hilbert term, Newton’s con¬ 
stant will become negative and gravity be a repulsive force. 
Certainly, this contradicts with our physical world, but this 
is not ruled out as a theoretical possibility. We know the 
Coulomb force is repulsive for like charges. 

The Hamiltonian of novel kinetic terms vanishes on-shell 
as well il, which is expected in covariant theories, so a 
bounded Hamiltonian is not a strong evidence for classical 
stability. One might need to examine other definitions of 
energy. As the local definition of gravitational energy is 
controversial, it may be more sensible to consider global 
energies (masses) according to the isometries of asymptotic 
spacetime. They are the conserved charges associated with 
the global symmetries of the vacuum where the infinite¬ 
dimensional diffeomorphism group is spontaneously broken 
to a finite-dimensional global symmetry group. In critical 
gravity models, the Abbott-Deser-Tekin mass of black 
hole solutions were shown to be zero EBl- As discussed 
before, the novel kinetic term £ 2 *^” can be considered as the 
special limit of critical gravity, and we expect it has the same 
property. 

From the quantum perspective, the correct sign for spin-2 
kinetic terms has a more definite answer. Unitarity requires 
particle poles to have positive residues. A propagator can be 
derived from the quadratic action, so the correct signs of the 
kinetic term and the mass term are determined. 
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The residue of a spin-2 propagator depends on the coupling 
to matter. It vanishes if the spin-2 field is not coupled to 
the energy-momentum tensor. So one can avoid negative 
residues by identifying the physical metric with the spin-2 
field with a correct sign kinetic term. Using the effective met¬ 
ric S, we can escape the problem of tree-level non-unitarity. 

Furthermore, to obtain the solutions of a model we need 
boundary conditions. When the Lagrangian allows for ghost¬ 
like excitations, they can still be avoided by proper bound¬ 
ary conditions. In this way, we could eliminate the ghost¬ 
like modes whose kinetic terms have wrong signs, then the 
bi-gravity models should reduce to healthy vector-tensor the¬ 
ories when the mass squared has a correct sign. 


VIII. CUT-OFF SCALE 


Ghosts are ubiquitous in the framework of effective field 
theories. Their presence do not disqualify the models from 
describing nature. They just tell us when the theories stop 
providing consistent descriptions and microscopic details be¬ 
come important. 

As an effective field theory, Einstein’s gravity has a cut-off 
scale set by the Planck mass. Higher curvature terms are also 
compatible with diffeomorphism invariance, so they should 
be present. The natural values of their coefficients are of or¬ 
der unity in terms of the Planck mass, and the cut-off scale 
remains the same. Ostrogradsky’s ghosts due to higher deriva¬ 
tive equations of motion are not excited below the Planck scale 
because their masses are around the Planck value. The cor¬ 
rections due to higher curvature terms are negligible at low 
energy scale. 

However, if for some unknown reasons, the coefficients of 
some correction terms are considerably larger than their nat¬ 
ural values, then we need to worry about Ostrogradsky’s in¬ 
stability even below the Planck scale. The cut-off scale of an 
effective field theory of gravity is lowered by the ghost modes. 

Let us consider an example in 4d which admits a bi-gravity 
reformulation 

E = lUp \/—g [^7? -t-77(e-t-ee-t-...) 

-|-A(1 -b e) -I- {ee + eee + .. (150) 

where the tensor structures are not written explicitly, the two 
spin-2 fields are dimensionless. After linearization 

and diagonalization, corresponds to the mass of the mas¬ 
sive spin-2 field which is ghost-like EE In principle, higher 
curvature terms are also allowed. 

Integrating out the auxiliary field e^i,, we have a La¬ 
grangian of higher curvature gravity, which schematically 


reads 


C = V^(^M^A 


Ml 


—Arr- 


(151) 


We can see a small mass in the bi-gravity formulation trans¬ 
lates into large higher curvature terms. 

To be more precise, there are two kinds of Ostrograd¬ 
sky’s ghosts in a model of higher curvature gravity, which 
can be rephased more transparently in the bi-gravity formu¬ 
lation. In the bi-gravity representation, the massive modes 
contain two kinds of ghost-like degrees of freedom. The first 
one is the ghost-like spin-2 mode due to a wrong sign kinetic 
term|3- The second one is Ostrogradsky’s scalar ghost or the 
Boulware-Deser ghost in a generic theory of massive spin-2 
field. 

Lor simplicity, let us assume Minkowski spacetime is the 
background solution. We also assume, after linearization and 
diagonalization, the kinetic terms of ( 1150b are given by the 
linearized Einstein-Hilbert terms to avoid more ghosts. These 
assumptions already constrain the possible terms in (1150b . 
The scale of the quadratic potential terms is set by the mass 
squared m^. Then both the spin-2 ghost and the scalar ghost 
are excited and interact with the healthy modes at a low scale 

A = m. (152) 


By requiring that the quadratic term of the linearized La¬ 
grangian takes the form of Lierz-Pauli mass term, the scalar 
ghost is absent in the quadratic action. But it can still appear 
in the interaction terms, which is known as the BD ghost. If 
we assume the effective Lagrangian is given by the bi-gravity 
models in the second class (l54l i. this scalar ghost can be elimi¬ 
nated completely. Then we can focus on the problem of spin-2 
ghosts. 

To increase the cut-off scale, let us first examine the 
quadratic action in detail. The massive spin-2 field is denoted 
by According to section|VI] its linearized action is given 
by Lierz-Pauli theory 

C = (153) 

where the kinetic term has a wrong sign and a is a model- 
dependent numerical factor. 

Let us decompose a la Helmholtz 


TT _ ^ TTl 


1 

Mp m 


{dpAy + d^Ap) + 


1 

Mp m? 


dpdi,cj}, 

(154) 


= df^Ap = 0. (155) 


It is assumed that the coefficients of the potential terms are of the same 
order, which is not a necessary assumption. When their magnitude are 
different, the smallest one is the most important. 


The kinetic term of the helicity-1 mode is from the mass terms, so they 
are ghost when the spin-2 ghost is also a tachyon. Let us assume the mass 
squared is positive. 
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where the dimensions of (p are 1 and the use of 

Mp, m is to canonically normalize the kinetic terms of the 
decomposed fields 1^ 

In terms of the decomposed modes, the Fierz-Pauli La- 
grangian becomes 


C=dp{H^),pd>^{H^r - 
+am2 


-2adpA, + 2a 


(158) 


where total derivative terms are neglected. The last term is a 
cross term, so we introduce 


Hp. = Hl, + ^<Pinp, (159) 

to diagonalize the kinetic terms. After diagonalization, the 
Lagrangian becomes 


The interaction terms start from the cubic order 

i + 2j + /e + Z = 3, 4, 5,..., (164) 

and a perturbative term 

((/.)'= (165) 

becomes important at the energy scale 

A = (166) 

The lowest scales of the interaction terms involving both the 
ghost-like modes (iT^, p) and the healthy modes (/i^, A) PI 
can be found: 

• The lowest scale of the two-derivative terms is 

A5 = (Mpm4)i/5, (167) 


C=dpH,pd^H''P - dpHpPd^H,'^ 

+am^ (5/77,"- 5^,5^") 

+ ‘^a^dp(j)d^(j) + 2c? rr?Hp^(j} + ^c?ir?(l? 

-2adpA^d>^A? (160) 

The kinetic terms of 5^, and (j> have wrong signs. If a > 0, 
then the helicity-1 modes Ap are healthy modes. Although 
the ghost-like modes can be excited at scale A = m, they 
are harmless before healthy degrees of freedom are coupled 
to them. The cut-off scales are then determined by the lowest 
scale of the interaction terms that involve both the ghosts and 
the healthy modes. 

By considering nonlinear redefinitions, we can always 
make the massless spin-2 field transverse 

1 rp 

hpv = T (161) 

Then we perform a general power counting of the possible 
perturbative terms without using the specific structures of the 
nonlinear terms. From the two derivative terms, we have 

(162) 

Assuming the coefficients of the potential terms in (I1501 l are 
of the same order, we have 

(163) 


where the cubic terms 

j = = 1, i = I = 0, (168) 

and 

k = 2, 1 = 1, i=j = 0 (169) 


become important. 

• The lowest scale of the potential terms is 


A3 = {Mpn?f/^, 

(170) 

where infinitely many terms 

j = l, fc = l-|-n, i = l = 0 

(171) 

and 

k = 2 + n, ( = 1, i = j = 0 

(172) 

II 

to 

(173) 


become important. 

The cut-off scale of the two-derivative terms is lower than 
that of the potential terms, so a generic bi-gravity model in the 
second class ( l54l i is a consistent effective field theory at least 
up to As. PI 


Another natural decomposition is with respect to the covariant derivative 
of the massless spin-2 field 






(156) 

(157) 


Mp m2 


= V('*)''Ap = 0. 


If a < 0, then the helicity-1 mode becomes a ghost due to a wrong sign 
kinetic term. If a = 0, the helicity-1 mode is strongly coupled. In addition, 
could be Ostrogradsky’s vector ghost if the equations of motion involve 
higher order time derivative terms of A^. But from the constraint analysis 
l3Qn . we can count the numbers of dynamical degrees of freedom, and we 
know the novel kinetic terms do not contain Ostrogradsky’s vector ghost. 
If the potential terms are modified, then Ostrogradsky’s scalar ghost is 
eliminated only at the linear level. If the massive spin-2 field has a correct 
sign kinetic term, then the cut-off scale set by the Boulware-Deser ghost is 
As as well. 
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We can further improve this by turning off then the 
kinetic terms have more gauge symmetries. The kinetic terms 
contain only helicity-2 modes whose the interaction scale is 
set by the Planck mass. In this way, we are able to increase 
the cut-off scale to at least A3 FI 

A ^ A3 = 1/3. (174) 

It should be noted that we need to make sure the spin-2 
ghost does not couple to the matter below the cut-off scale. 
This indicates we should consider an effective metric ll^ FI 

In a different region of parameter space, it is possible that 
the massless spin-2 field has a wrong sign, while the massive 
one is healthy. By eliminating the BD-ghost and using the 
gauge-invariant kinetic terms, the cut-off scale is set by the 
interaction terms k = 2 + n, I = 1, i = j = 0 from the 
potentials, which is A3 again. When the mass squared has a 
correct sign, these effective field theory contain a ghost-free 
massive graviton and a decoupled ghost-like, massless spin-2 
field below the cut-off scale. 

Here we want to give one possible reason for “naturally" 
large coefficients of the higher curvature terms. In the bi¬ 
gravity formulation, the gauge symmetries are enhanced when 

= 0, 03 = 0, (175) 

so small values of m? and 03 are technically natural, which 
is analogous to the mass of electron. Quantum corrections to 
these parameters should be multiplicative, rather than addi¬ 
tive. 


IX. SYMMETRIC CONDITION 


(179) 

(180) 

(181) 

We argue that the equations of motion for generally 

lead to the symmetric condition 

because the equations for can be written in a matrix form 

=0, (183) 


which gives the symmetric condition if A is invertible. An 
important point is that the Lagrangians do not contain linear 
terms, so the equations of motion for are homogeneousEB 
The argument is clear if (11811) is not considered. When the 
Lagrangian contains (11811) . the equations of motion contain 
cubic terms of Then we can write the cubic terms as 
products of quadratic terms and linear terms, and think of the 
quadratic terms as part of the matrix A. El 

We do not have a proof that the above argument works in 
general, but we check several examples and always find that 

det.A7^0. (184) 

The spirit is close to Hi , where the symmetric condition is 
derived from a local Lorentz transformation. In addition, we 
do not rule out the possibility that A could be degenerate at 
some singular points of the phase space. 


In the above sections, we impose the symmetric condition 
or the Deser-van Nieuwenhuizen condition S to minimalize 
the numbers of dynamical degrees of freedom. In this section, 
we present a general way to derive the symmetric condition 
from the equations of motion. 

Let us decompose the rank-2 tensor e^i, into two parts 

ep. = F^^E.^riAB = (176) 

where is the symmetric part and is the antisymmetric 
part 




= e 


I'll ’ 


'-' 111 ' 


= —e 


Ufl ■ 


(177) 


In the Lagrangians, an antisymmetric product vanishes if it 
contains an odd number of In 4d bi-gravity models, the 
antisymmetric part of ep,y only appears in the terms below 


(178) 


3'* In Isil) . it was shown that the A 3 interaction terms vanish in some cases 
which are equivalent to quadratic curvature gravity. This indicates the cut¬ 
off scale may be higher than A 3 . 

But the BD ghost will be revived at some scale above A 3 due to the use of 
an effective metric. 


X. HIGHER DERIVATIVE GENERALIZATIONS 


Along the line of Lovelock terms, the novel kinetic terms 
can be generalized to novel higher-derivative terms 101 El 

R{E) A • • • A R{E) AEA---AEAEA---AE, (185) 

which might be inconsistent with terms involving R{F). 

In section rvl we discuss the special cases with only one F 
vielbein 

R{E) A---AR{E) AE A---AE AF, (186) 


We assume the matter does not couple to the antisymmetric part of epi, 
linearly. For example, if the physical vielbein is a linear combination of 
and , the corresponding physical metric will contain a quadratic 
term so the equations of motion for are homogeneous. 

Note that A can be degenerate for special values of if they are not the 
solutions of the equations of motion at the same time. 

If we impose the symmetric condition and fix the second metric to be 
Minkowski, they are equivalent to the higher derivative interactions pro¬ 
posed in (H. 
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which have the same gauge symmetries as R{E) A E A F. 
They describe the derivative interactions between two mass¬ 
less, gauge invariant spin-2 fields, where the Boulware-Deser 
ghost is absent. We expect other higher derivative terms in 
(11851) do not contain the BD ghost as well. 

For the extension to multi-gravity, we have 

R{E^^^) A • • • A R{E^^'>) A E^^'> A • • • A (187) 

where n is the number of curvature 2-forms and E^^'> vieT 
beins can be the same or different. Lovelock terms and dRGT 
terms are unified in (11871) . 

In the end, we would like to connect with some results in 
the literature. The bi-gravity models with (11861) in the met¬ 
ric formulation were already proposed in Oil] as generaliza¬ 
tions of New Massive Gravity. The BD ghost was argued to 
be absent by counting the degrees of freedom using symme¬ 
tries OHd^ . The antisymmetric structure guarantees that the 
equations of motion for the decomposed fields 

-f- S/ 4- SI -f SI (188) 

are of second order, so additional degrees of freedom are 
avoidedl^. Then one can count the dynamical degrees of free¬ 
dom in the bi-gravity models and show that the total number 
is at most {<P — 2d — 1), which is that of a massless and a 
massive spin-2 fields. Therefore, the Boulware-Deser ghost is 

absent0 

From this argument, we can see why the curvature tensors 
should be associated with the same spin-2 field If a 

curvature tensor contains the second spin-2 field then the 
equations of motion for the decomposed modes will usually 
be of higher order, because they are not gauge modes in a 
curvature tensor and no apparent antisymmetric structure is 
protecting them EH- 

The decomposed field argument concerning ( 11861 ) is based 
on the fact that Lovelock tensors are divergence-free. For the 
other novel derivative terms in (11871) . we can generalize this 
argument by using the second Bianchi identity 

= (189) 

The covariant derivatives in front of the decomposed fields 
will not act on the Riemann tensor after integrating by parts. 


This ai'gument is dangerous. The equation of motion for the decomposed 
field (p contains third order derivative terms of the metric, in the form of co¬ 
variant derivatives of curvature tensors. In addition, if one varies the action 
with respect to the metric after the substitution, the equations of motion will 
contain third order derivative terms of (p due to the variations of covariant 
derivatives. But it is possible that the third order time derivative terms can 
be removed by the time derivatives of some second order equations f^ . 
then the counting of the degrees of freedom is correct. 

When there are additional gauge symmetries, the number of dynamical 
degrees of freedom is reduced. 

In (54h . an Einstein-Hilbert term for was introduced to obtain unitary 
models. However, we suspect the absence of ghost-like degrees of freedom 
is an artifact of linearization. For example, in the minisuperspace approxi¬ 
mation, the Hamiltonians are not linear in the lapse functions. 


SO the equations of motion for the decomposed fields will not 
contain 4th order derivative terms of the metric. The variation 
of a Riemann tensor contains some second covariant 

derivative terms of which are antisymmetrized , so the 
equations of motion for the metric will not contain 4th order 
derivative terms of (j). 

In the vielbein formulation, the second Bianchi identity 
stems from a basic identity of exterior derivative 

d? = 0, (190) 

which is the key element of the unifying framework 10,[Ml- 


XI. CONCLUSION 

In summary, we present evidence that 

RiE^^P A • • • A R{E^'^p A E^‘^'> A • • • A (191) 

are basic building blocks for the actions of interacting spin- 
2 fields that are free of the Boulware-Deser ghost 0 Mod¬ 
els that can be constructed from these building blocks include 
Einstein gravity, Weyl gravity. Lovelock gravity. New Mas¬ 
sive Gravity, dRGT massive gravity and some of their gener¬ 
alizations. The parameter space is further extended by novel 
derivative terms. 

Curiously, the building blocks (11911 can be obtained from 
Lovelock terms by replacing some of the vielbeins in the 
wedge products with other vielbeins. 

The novel two-derivative terms in 4d are studied in detail: 

• Based on a minisuperspace analysis, a large class of bi¬ 
gravity models (l54l) are identified, which are potentially 
free of the Boulware-Deser ghost. 

• The bi-gravity models in this class (l54l) do not have the 
usual single dynamical metric limit with a fixed metric, 
which is in accordance with the no-go theorem for new 
kinetic interaction for single dynamical metric in fTM . 

• We reformulate some well-understood models of higher 
curvature gravity as bi-gravity models in this class ( |5^ , 
Their spectra are known to contain 1 massless and 1 
massive spin-2 fields without the BD ghost. 

• The argument that New Massive Gravity is free of the 
BD ghost is extended to other bi-gravity models in 
this class (l54]). w hich applies to novel higher-derivative 
terms as wellrd. 


Parity is assumed to be preserved, otherwise there are more possible terms. 
For example, in 3d, one could introduce a gravitational Chem-Simons term 
which violates parity (5^. The critical points of higher derivative gravity 
theories were first investigated in this context l5^ . 

The absence of the B D g host in 4d novel kinetic terms is proved hy the 
constraint analyses in l30ll . 
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This class of bi-gravity models are interesting despite the 
issue of spin-2 ghosts. Firstly, as toy models of quantum grav¬ 
ity, they have better chance to be perturbative renormalizable 
and there are less negative norm states because the BD-ghost 
is absent. Secondly, as effective field theories of gravity, they 
can increase the cut-off scale set by higher derivative terms 
with large coefficients PI 

In general, we can avoid the ghost modes by reducing the 
number of dynamical degrees of freedom. A useful strategy 
of eliminating ghost modes is to impose specific boundary 
conditions. Another method to remove the spin-2 ghosts is 
simply setting the decomposed helicity-2 modes in (1188b 
to zero. They may give rise to healthy vector-tensor theories 
0 

In the light of the AdS/CFT correspondence ifs^ . the large 
class of non-unitary bi-gravity models may be dual to non¬ 
unitary conformal field theories. It is interesting to explore 
non-unitary holography in the extended parameter space. 
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